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We study the dynamics of strongly repulsive Bose gas in tilted or driven bichromatic optical 
lattices. Using the Bose-Fermi mapping and exact numerical method, we calculate the reduced 
single-particle density matrices, and study the dynamics of density profile, momentum distribution 
and condensate fraction. We show the oscillating and breathing mode of dynamics, and depletion of 
condensate for short time dynamics. For long time dynamics, we clearly show the reconstruction of 
system at integer multiples of Bloch-Zener time. We also show how to achieve clear Bloch oscillation 
and Landau- Zener tunnelling for many-particle systems. 
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I. INTRODUCTION 

The dynamics of a particle in a period structure has 
been a fundamental subject, with the eigenenergies form- 
ing the famous Bloch bands [1] and eigenstates being 
delocalized. If a weak external constant force is intro- 
duced, contrary to our intuition, the particle undergoes 
an oscillatory motion rather than uniform motion due 
to acceleration by the force, which is known as the fa- 
mous Bloch oscillation [l|, Q . Under single-band tight- 
binding approximation, eigenenergies of the system form 
the Wannier-Stark ladder [3] with eigenstates localized. 
Bloch oscillation bas been observed in semiconductor su- 

t)erlattices for electrons Q , optical lattices for cold atoms 
2] and photonic crystals for light pluses [El, Hf. For a 
stronger force, a directed motion is re-introduced by re- 
peated Landau-Zener tunneling to higher Bloch bands 
[THIGI. For usual cosine-shaped potentials, band gaps 
decrease rapidly as the energy increases, and this would 
lead to the decay of Bloch oscillation, which has been 
observed in Ref . [iMl . 

In order to study the interplay between Bloch oscilla- 
tion and Landau-Zener tunnehng, one needs at least a 
two-band system with the lowest two bands being well 
separated from upper ones. Furthermore, the gap be- 
tween the lowest two Bloch bands must be small for ob- 
serving clear signal of Landau-Zener tunneling. This can 
be achieved by bichromatic optical lattices [14], where 
parameters of the system are adjustable and control- 
lable. Bichromatic lattices have been implemented by 
superimposing two incoherent optical lattices with the 
wavelength of one lattice being two times of the other 
[l6| , o r b y virtual two-photon and four-photon pro- 
cesses [l7l-[21|. Under two-band tight-binding approxi- 
mation, eigenenergies form two Wannier-Stark ladders 
with energy spacing doubled and an offset between them 
[22|,[23j. The corresponding eigenstates are still localized. 
The dynamics of a particle is governed by two timescales, 
i.e., the Bloch period and the period of Zener oscillation. 
If the two periods are commensurate, system reconstructs 
at integer multiples of Bloch-Zener time. 



So far, most of works concentrate on the dynam- 
ics of the single-particle system. The generalization of 
these results to interacting many-particle systems re- 
mains an open question. The possibility to investigate 
Bloch oscillation and Landau-Zener tunneling of inter- 
acting Bose-Einstein condensate (BEC) experimentally 
has attracted much interest. Most of theoretical studies 
are based on the mean-field approximation (241-1271] and 
Gross-Pitaevskii (GP) equations [28]. Results for the dy- 
namics of strongly interacting many-particle systems are 
rarely known. In this paper, we study the dynamics of 
interacting bosons in bichromatic optical lattices under 
constant drag force in the limiting case with infinitely 
repulsive interaction which permits us to solve the prob- 
lem exactly. The one-dimensional (ID) Bose gas with 
infinitely repulsive interaction is known as the hard-core 
boson (HCB) or Tonks- Girardeau (TG) gas [29], which 
can be exactly solved via the Bose-Fermi mapping l29l 
and has attracted intensive theoretical attention |3Q|-132| . 
Experimental access to the required parameter regime 
has made the TG gas a physical reality (33|, [s^ • Follow- 
ing the exact numerical approach proposed by Rigol and 
Muramatsu [35, 36] , we calculate the dynamics of density 
profile, momentum distribution, and condensate fraction 
for hard-core bosons in the tilted bichromatic optical lat- 
tice, and show Bloch oscillation, Landau-Zener tunneling 
and reconstruction of the system at integer multiples of 
Bloch-Zener time. 



The paper is organized as follows. In Section II, we 
present the model and the exact approach used in this 
paper. We also recover the dynamics of single parti- 
cle system in this section. In Section III, we study the 
short-time and long-time dynamics for hard-core bosons 
in bichromatic optical lattice with a constant drag force. 
We also show how to achieve clear Bloch oscillation and 
Landau-Zener tunneling for the many-particle system. 
Finally, a summary is presented in Section IV. 
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II. MODEL AND METHOD 

In the present section we describe the exact approach 
which we used to study ID hard-core bosons in tilted 
or driven bichromatic optical lattices. Under the tight- 
binding approximation, the system can be described by 
the following Hamiltonian: 

H = - J^(6|6,+i + H.c.) + S Y,{-iyn^i 



(1) 



Here we only consider the nearest-neighbor hopping and 
neglect the off-diagonal terms of position operator x in 
the Wannier basis. The operator {hi) is the creation 
(annihilation) operator of boson which fulfills the hard- 
core constraints [35], i.e., the on-site anticommutation 
h\} = 1) and [6^, h^-] = for z 7^ j; n\ is the bosonic 
particle number operator; J is the hopping amplitude 
being set to be unit of energy (J = 1); Vi = Vnii — io)'^ 
is the harmonic potential for preparing the initial state of 
system, with Vh the strength and io the position of the 
trap center; S is the energy shift of alternate site; F is 
the strength of driven force. For convenience, the lattice 
spacing is set to be unit. 

In order to study the dynamics of the hard-core bosons 
in driven optical lattice, we first load the hard-core 
bosons into a bichromatic optical lattice with an addi- 
tional harmonic trap, then we switch off the harmonic 
trap and turn on the driven force. We shall study the 
evolution of the initially prepared system and dynamics 
of the system under the driven force. First, the initial 
state is the ground state of Hamiltonian: 

Hinit = -jY,{b\hi+i + H.C.) + ,5^(-l)X^ + 

i i i 

(2) 

with particle number N. In order to get the initial state, 
it is convenient to use the Jordan- Wigner transformation 
[13 (JWT): 



/in< 



/3=1 



(3) 



to map the Hamiltonian of hard-core bosons into the 
Hamiltonian of noninteracting spinless fermions H^:^^^ 
which is in the same form as i^init •, but with all the boson 
operators, e.g.,6j, hi and n^, being replaced by the cor- 
responding fermion operators, e.g., //, fi and n{ . The 
ground-state wave function of the system with N spinless 
free fermions, which is a product of lowest N eigenfunc- 
tions, can be obtained by diagonalizing H^:^^ and can be 
represented as 



N 



TT^Pi„//|o), 



(4) 



where L is the number of lattice sites, N is the number 
of fermions (same as bosons), and coefficients Pin are the 
amplitude of the n-th single-particle eigenfunction at z-th 
site which can form an L x matrix P [ss]. 

After releasing from the trap, the system is described 
by Hamiltonian: 



= -J J2iblbi+i + H.c.) + 6 yj-iyn'^ 



i i i 

(5) 

Similar to the above method, from the corresponding 
free- fermion Hamiltonian , we can get all single- 
particle states and corresponding energies, and use 
to represent all the single-particle states, Si to represent 
the energies. The nonequilibrium quantum dynamical 
properties of system can be calculated through the equal 
time one-particle Green function which is defined as: 



hcb(*)), 



(6) 



with I^hcb(^)) is the wave function of hard-core bosons 
at time t after releasing from harmonic trap. After some 
derivations (see Appendix A), one can get 



Gij{t)=det \{P^y P^ 



(7) 



where P^ and P^ are obtained from P, P' and Si. It 
follows that the reduced single-particle density matrix 
can be determined by the expression 

Pij{t) = {h\hj)t = Gji(t) + 5ij(l - 2Gii{t)). (8) 

The momentum distribution is defined by Fourier trans- 
form with respect to i — j of the reduced single-particle 
density matrix with the form 



-ik{i-j) 



(9) 



where k denotes momentum. The natural orbitals (j)^ are 
defined as eigenfunctions of the reduced single-particle 
density matrix (soj . 



(10) 



=1 i=i 



The natural orbitals can be understood as being effective 
single-particle states with occupations A^. For noninter- 
acting bosons, all particles occupy in the lowest natural 
orbital and bosons are in the BEG phase at zero temper- 
ature, whereas only the quasicondensation exists for ID 
hard-core bosons [35] . 

For hard-core bosons we know that the state for the 
corresponding Fermi Hamiltonian is a product of time- 
dependent single-particle states and each single-particle 
state evolves itself (see Eq. ()A4p in Appendix A). We re- 
cover the single-particle properties of the Hamiltonian 
He in Appendix B, which have been studied by Breid 
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et. al [22|. From Appendix B, we know that a single- 
particle wave function will be reconstructed at integer 
multiples of Bloch-Zener time (Tbz), if Ti and T2 are 
commensurate. Here Ti is the Bloch period decided by 
the spacing of Wannier- Stark ladders, and T2 is the pe- 
riod of Zener oscillation decided by the offset between 
Wannier- Stark ladders. Since each single-particle state 
can be reconstructed at period of time, and the period 
is independent of the single-particle state and decided 
by parameters of Hamiltonian, the state for many-body 
hard-core bosons composed of a product of single-particle 
states will also be reconstructed at integer multiples of 
Bloch-Zener time. 



III. QUANTUM DYNAMICS OF HARD-CORE 
BOSONS 

In order to observe reconstruction of system, two peri- 
ods Ti and T2 must be commensurate which are decided 
by F and where £^0 is a function of F and 5. In 
FigJTJi, we show numerical results of £^0 versus 5 for a 
particular F. For different F, structure of the picture is 
similar. In order to generate particular Bloch-Zener time 
(Tbz), ^ must be one of the discrete numbers. For exam- 
ple, if we want Tbz = Tb for the system with F = 0.05, 
we have to let T^o = and then 6 = 0.1684, 0.3614, .... 

For comparison, we first recover the dynamics of a sin- 
gle particle system. We choose Eq = and let system 
reconstruct at integer multiples of Bloch time {Tb). From 
now on, we take Bloch time Tb = 2Ti as the reference 
timescale. In Fig(Tl3, we show the dynamical evolution 
of density profile for a single-particle system, from which 
one can see the reconstruction of the density profile at 
integer multiples of Bloch time. The edge of Brillouin 
zone is reached at t = Tb/^^ and part of particle moves 
into upper excited Bloch band located in upper half of 
the figure. The particle in upper excited band returns to 
lower Bloch band at t = 3Tb/4. For quantitative analy- 
sis of Landau-Zener tunneling rate, one can characterize 
it by the number of particles in upper half of picture at 
time t = Tb/'^. Numerical results show that the Landau- 
Zener tunnelling probability [2 81] : 

PLZ«exp(-^). (11) 

So in order to see a clear signal of Landau-Zener tun- 
neling, we have to choose small 5 for a given strength of 
force F. Furthermore, the available interval for motion 
of a particle in lattice [^] : 

L ^ 4./F. (12) 

In FigHt, we also show the dynamical evolution of mo- 
mentum distribution of the single-particle system. Par- 
ticles with momentum in the interval (— 7r/2,7r/2) are in 
lower Bloch band and outside the region particles are in 
upper excited Bloch band. From this picture we can see 
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FIG. 1: (Color online) (a): ^0 vs 6 for the system with F — 
0.05. The dynamics of density profile (b) and momentum 
distribution(c) for single particle system with 5 — 0.1684, F — 
0.05, and Vh — 0.001. (d): The breathing mode dynamics 
of density profile for single particle system with 5 = 0.1684, 
F = 0.05, and Vh = 20. 



the clear signals for Bloch oscillation and Landau-Zener 
tunneling. Furthermore the momentum is linear with 
time with slope given by F. Next we consider a localized 
initial state which has a wide momentum distribution and 
can be implemented by setting a strong harmonic trap 
potential, for example, Vh = 20 here. From Fig(l]i, one 
can find a breathing behavior of density profile. The re- 
construction still happens at integer multiples of Tb. But 
one can observe that the enveloping structure with period 
of Tb is overlayed by a breathing mode of smaller am- 
plitude, whereas part of particle remains in lower Bloch 
band all the time with period of Te/2. 

Now we consider the case of many-body hard-core 
bosons and study the short time dynamics firstly. Sys- 
tems with various particle numbers will be considered 
by keeping the other parameters fixed. For comparison, 
we show the dynamics of a single-particle system in the 
first column of Figl2j After releasing from trap, the par- 
ticle speeds up under the drag force T, and it reaches 
the edge of Brillouin zone at time Tb/4. Part of particle 
moves into upper excited Bloch band through Landau- 
Zener tunneling, while the other part of particle remains 
in lower Bloch band with changing of sign of momentum 
by Brag scattering. Because of Landau-Zener tunneling, 
the particle turns into two parts and they are separated in 
the real space with particles in the upper half of the den- 
sity distribution being in the upper excited Bloch band. 
In momentum distribution, the particles outside the first 
Brillouin zone of (— 7r/2,7r/2) are in upper excited Bloch 
band. The changing of sign of momentum occurs at time 
Tb /4 by Brag scattering for particles in the lower Bloch 
band. In the third row of Figl2l we show the evolution of 
condensate fraction which is defined as Aq/A^ with Aq be- 
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N=l N=2 N=5 N=30 N=150 




FIG. 2: (Color online) The short-time dynamical evolution of density profile (first row), momentum distribution (second row), 
condensate fraction (third row), and modules of the reduced single-particle density matrix (fourth row) at time t = Tb/2 for 
systems of different particle number with F = 0.05, 6 = 0.1684, Eq = 0, and Vh = 10~^ (Vh = 8 x 10~^ for the fifth column). 



ing the occupation of the lowest effective single-particle 
state. For single particle system, the particle is always 
in the lowest effective single-particle state, and the con- 
densate fraction is one all the time. For the ID many- 
body systems of hard-core bosons, there are only quasi- 
condensation with Aq cx a/]V [29|, [ssj . In the fourth row of 
Figl2l we show the reduced single-particle density matrix 
at time Tb/2. Here we consider the modulus because of 
the elements of density matrix being complex numbers 
after turning off trap. The upper right spot in picture is 
caused by particles in upper excited Bloch band and the 
lower left spot for lower Bloch band. For single particle 
dynamics, however particles in upper excited Bloch band 
and lower Bloch band are separated in real space, there 
are phase coherence between them, the off-diagonal parts 
of reduced single-particle density matrix are very strong. 

In the second to fifth column, we show the dynamics 
for the hard-core bosons with N = 2,5,30 and 150, re- 
spectively. As particle number increases, the adding par- 
ticle has to occupy higher single-particle state because 
one state can only be occupied by a hard-core boson. 
And the size {S) of system becomes larger and larger, 
while the available interval L of the system decided by 
force F remains unchanged. Also, the width of momen- 
tum distribution becomes wider for the larger system. 
As the momentum distribution becomes wider, it takes 
shorter time for particles at the edge of the momentum 
distribution to reach the edge of Brillouin zone, and thus 
Landau- Zener tunneling happens early, which leads to 



the condensate fraction decreasing early. The conden- 
sate fraction decreases (cx I/a/TV) as the particle number 
increases. As time increases but is smaller than Te/4, 
the condensate fraction basically does not change. A 
slight increase of the condensate fraction in short time 
is caused by the expanding after turning off trap [38j. 
At time t = Tb/^^ Landau- Zener tunneling happens, 
and the condensate fraction decreases quickly. After this 
there appears an overdamped area, and then the conden- 
sate fraction keeps unchanged basically. For the reduced 
single-particle density matrix at time Te/2, lengths for 
the two parts of diagonal terms become larger as the size 
(yS) of system increases. The off-diagonal parts of matrix 
become weaker as particles add in. Also, the phase coher- 
ence between particles in the upper excited Bloch band 
and lower Bloch band decreases. As shown in the figure, 
the particles between upper band and lower band lost 
their phase coherence when = 30. Furthermore, the 
particles in the upper excited Bloch band and lower Bloch 
band developed phase coherence inside each part, the re- 
duced single-particle density matrix has exponential-law 
decay in each part as the distance increases. 

We note that the density profile no longer splits into 
two obviously separated parts after time Te/4 for sys- 
tem with = 30 as shown in the fourth column, where 
the particles in the upper excited Bloch band and lower 
Bloch band are overlapped in real space. To achieve this 
situation, one has to adjust the parameters of system 
and let the size of system larger than the available in- 



terval {S > L). Also one has to avoid the initial state 
staying in a localized state because localized initial state 
will cause the breathing-mode dynamics. Although the 
two parts are overlapped in real space, they are sepa- 
rated in the momentum distribution. The width of the 
momentum distribution becomes larger as the particle 
number increases. The dynamical evolution of conden- 
sate fraction still has similar structure, except that the 
quick decrease happens more early. As two parts of par- 
ticles are overlapped in real space, the diagonal terms of 
density matrix are also overlapped. For sites outside the 
overlapped area in upper excited Bloch band or lower 
one, the density matrix still has a exponential-law de- 
cay as distance increases. For overlapped area, the den- 
sity matrix is irregular, but as the distance increases, 
it goes to zero quickly. In the fifth column, we show 
the dynamics of system with = 150. As the particle 
number increases, the localization of initial state becomes 
stronger [35], and we see the breathing mode of dynam- 
ics of density distribution. Lots of particles are localized 
in the center, and particles in two parts are overlapped. 
The dynamics of momentum distribution is also changed. 
The momentum distribution does not increase linearly as 
time increases, and momentum distributions for particles 
in the upper excited Bloch band and lower Bloch band 
are overlapped. For localized initial state, the momen- 
tum distribution is almost flat. Right after turning off 
trap, there are particles moving into upper excited Bloch 
band through Landau-Zener tunnelling, and the conden- 
sate fraction decreases immediately. Meanwhile, for a 
localized initial state, the expanding is also important 
after turning off trap, and there is a large over damped 
area in the breathing mode dynamics of condensate frac- 
tion. The reduced single-particle density matrix still has 
exponential-law decay for particles in localized states. 

Next, we consider the long-time dynamics of hard-core 
bosons. In Figl3^, we show the dynamical evolution of 
density distribution for five-particle system with Bloch- 
Zener time Tbz = Tb, where Eq = 0. After releasing 
from the harmonic trap, particles move along with the 
direction of the force F. Around time t = Tb/^^ par- 
ticles reach the edge of Brillouin zone, and part of par- 
ticles move into the upper excited Bloch band through 
Landau-Zener tunneling and keep moving along with 
force F. The other part of particles remain in the lower 
Bloch band, but the momentum of particles changes the 
sign due to Brag scattering and thus the particles move 
against the force F. Around time t = 3Te/4, particles 
reach the edge of Brillouin zone again, and all particles 
move into the lower Bloch band. At time t = T^, the 
system reconstructs into the initial state. As time goes 
on, more periods occur. The dynamical evolution of the 
momentum distribution for the same system is shown in 
FiglSb. After turning off trap, particles accelerate under 
force F, and the momentum increases linearly with time. 
Around time t = Tb/^^ particles reach the edge of Bril- 
louin zone {k = 7r/2). Part of particles, which move into 
the upper Bloch band through Landau-Zener tunnelling. 
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FIG. 3: (Color online) The dynamical evolution of density 
profile (a,d), momentum distribution (b,e), and condensate 
fraction (c,f) for five-particle system with F — 0.05, Vh — 
10"^ and 5 = 0.1684 (a,b,c), 5 = 0.238 (d,e,f). 



remain in the same belt and their momentum increases 
linearly with time. The other part of particles keep in 
the lower Bloch band and change the sign of their mo- 
mentum by Brag scattering. Two parts of particles are 
recombined around time t = 3Te/4, and the momentum 
distribution returns to the initial distribution at t = Tb- 
In FiglSt, we show the dynamical evolution of conden- 
sate fraction. Around time t = Tb/^^ the condensate 
fraction has a quick decrease due to Landau-Zener tun- 
neling. After this there is an overdamped area. Around 
time t = 3Te/4, two parts of particles recombine, and 
the condensate fraction increases quickly. The system 
returns to its initial state at time T = Tbz- Further- 
more, the condensate fraction is symmetrical with center 
at T = Tbz/2. In Fig 131 we also show the dynamics of 
system with Tbz = 2Te where Eq = T/4, which has sim- 
ilar properties with the previous one. From Figl3]i, one 
can obviously observe the reconstruction of system af- 
ter Bloch-Zener time Tbz, Bloch oscillation and Landau- 
Zener tunneling. 

In FiglH we show the dynamical evolution of sys- 
tem with eighty hard-core bosons. As particles add in, 
the particles in the initial state get more localized, and 
the dynamical evolution of the system is dominated by 
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FIG. 4: (Color online) The dynamical evolution of density 
profile (a,d), momentum distribution (b,e), and condensate 
fraction (c,f) for the eighty-particle system with F = 0.05, 
S = 0.1684 and Vh = 10"^ (a,b,c), Vh = 4 x 10"^ (d,e,f). 



the breathing mode instead of the oscillating mode. In 
FigHH, we show the dynamical evolution of the den- 
sity profile for system with F — 0.05, (5 = 0.1684 and 
Vh — 10~^. This picture is not in prefect breathing mode 
because of the initial state being not localized enough. 
First of all, the reconstruction of density profile happens 
again after integer multiples of Bloch-Zener time Tbz- 
Second, lots of particles are localized at the center and 
form a belt. Third, one can observe that the enveloping 
structure is overlayed by a breathing mode of smaller am- 
plitude. The outside breathing mode is formed by Bloch 
oscillation and Landau-Zener tunneling for particles in 
the upper excited Bloch band with period Tbz = T^b- 
The inside breathing mode is formed by Bloch oscillation 
for particles remaining in the lower Bloch band with pe- 
riod Ti = Tb/2. Center is the localized belt. In FigI3]3 we 
show the dynamical evolution of the momentum distribu- 
tion. Strip structure disappears and windmill structure 
appears with the center at (nT^, /c = 0)(n G N). In this 
picture we can not distinguish the two parts of particles, 
and lots of particles are localized at area of A: ^ all the 
time. Foremost, the momentum distribution reconstructs 
at integer multiples of Tbz- The dynamical evolution of 
the condensate fraction is shown in FigHt, As the initial 





FIG. 5: (Color online) The dynamical evolution of density 
profile (a), momentum distribution (b), for the eighty-particle 
system with F 0.01, b 0.08211 and = 4 x 10~^ 



state is a localized state, there are many particles with 
high momentum, and right after turning off trap the con- 
densate fraction decreases quickly, and it reconstructs at 
t — Tbz- As a lot of particles add in, the initial state 
becomes a localized state and one can observe breathing 
mode dynamical evolution of density profile. In order to 
observe the oscillating mode and Landau-Zener tunnel- 
ing, we have to reduce the localization of initial state by 
decreasing the strength of the harmonic trap. The dy- 
namics of system after decreasing the strength of trap is 
shown in the second column of FigHJ For the density 
distribution we actually see the oscillating mode, but we 
can not see the Landau-Zener tunneling clearly. At time 
range (Tb/4, 3Tb/4), the particles in the upper excited 
Bloch band and lower Bloch band are overlapped in real 
space. This is due to the available interval for the motion 
being L ^ 4/F = 80 for system with F = 0.05, but the 
size of system is about S ^ 300. It is clear that S > L 
and the overlapped structure appears. After decreasing 
the strength of trap, the strip structure reappears in the 
dynamical evolution of the momentum distribution in- 
stead of the windmill structure. Furthermore the mo- 
mentum distribution between the particles in the upper 
excited and lower Bloch band are distinguishable despite 
they are overlapped in real space. For the condensate 
fraction, the curve is fiat again for the time short after 
turning off trap. 

So in order to observe the oscillating mode of dynamics 
and Landau-Zener tunneling, one has to let L > S. Now 
we have to decrease the strength of force F to achieve 
a bigger available interval for the motion. Once F is 
changed , we have to change 6 too to make £^o = 
for the system still with Bloch-Zener time Tbz = Te- 
Furthermore we have to choose a small S to achieve the 
big enough Landau-Zener tunneling probability, other- 
wise we only see the Bloch oscillation but can not see 
the Landau-Zener tunneling. In FiglS^, we show the dy- 
namical evolution of density profile for the system with 
F = 0.01, S = 0.08211 and Tbz = T^. Now we see the 
clear signals of Bloch oscillation and Landau-Zener tun- 
neling. The dynamical evolution of the momentum dis- 
tribution and condensate fraction are similar to Fig|4^ 
and FigHF. 
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IV. CONCLUSION 

In summary, we have studied the dynamics of infinitely 
repulsive Bose gas in tilted or driven bichromatic opti- 
cal lattices. Using the Bose- Fermi mapping and exact 
numerical method, we calculate the one-particle density 
matrices, density profiles, momentum distributions, nat- 
ural orbitals and their occupations (condensate fraction). 
Both the short-time and long-time dynamical evolution 
of density profile, momentum distribution and conden- 
sate fraction are studied. It is clearly shown that the 
reconstruction of system at integer multiples of Landau- 
Zener time. We also give estimations for how to achieve 
clear Bloch oscillation and Landau- Zener tunneling in 
given many-particle systems. 
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Appendix A: THE EQUAL TIME GREEN 
FUNCTION FOR HARD-CORE BOSONS 



The equal time Green function for the hard-core bosons 
can be written in the form 



G^J{t) = (^HCB(t)|6,6]|^HCB(t)) 

(3=1 7=1 

= (^^1^^) (Al) 

where |^hcb(^)) is the wave function of hard-core bosons 
at time t after releasing from harmonic trap and 
is the corresponding one for noninter acting fermions. In 
addition, we denote 

i-i 



7=1 
i-1 



/3=1 



(A2) 



The wave function \'^F{t)) can be easily calculated 
with the initial wave function |^^) 



N L 



= e-^^= *|*^) = n E Pin{t)fl\0), (A3) 



n=l 1=1 



with 



N N 

PUt) =Y.^-''''PlkY.^PTkPjn). (A4) 

fe=l 3=1 



where we have set = 1 in evolution operator, and P{t) 
is the matrix of \^ pif)) in the same way as |^^). In 
order to get Eq. ()A3[) . one has to insert XljLi I — ^ 

in to it, where = Ylt=i -^nj fn\^) lowest j-th 

eigenfunction of .We can see that is still a 

product of time-dependent single-particle states. 
In order to calculate ^^(and we notice that 



i-1 

7=1 



J-1 



(A5) 



7=1 



Then, the action of Y^^^i^^'"''^^^^^ on the state |^f(^)) 
(Eq. ()A3p ) generates only a change of sign on the element 
Pin{t) for I < and one has to add a column to P{t) 
with element Pj,N-\-i = 1 and all the others equal to zero 
for the further creation of a particle at site j. Then 



N+l L 

n=l ^=1 
AT+l L 

n=l 1=1 



(A6) 



where and P^ are obtained from P{t) changing the 
required signs and adding the new column. 
The Green function is written as 



AT+l L 



N+l L 



n=l 1=1 
L 



n' = l I' 



^hl ' ' ' ^In+iN+1 

/l---l7V + l,^l---l'7V + l 

= det[{P^)^P^] (A7) 
which requires 

(o|/h-- fi.+Jl ■■■f}'\o) = 

'■N + l ^ 

(A8) 

with the Levi-Civita symbol and A = l-'-N + l. 



Appendix B: THE SINGLE-PARTICLE 
PROPERTIES OF HAMILTONIAN He 

First of all, for the field- free case with F = 0, a 
straightforward calculation yields the dispersion relation 



(3k 



(-l)^+V^' + 4cos2(fe), 



(Bl) 



and corresponding wavefunctions \x(3{k)) (Bloch bands 
and Bloch waves) with the miniband index /3 = 0, 1. For 
nonzero F, the spectrum of Hamiltonian consists of two 



Wannier- Stark ladders with an offset in between. After 
introducing translation operator 



(B2) 



n— — oo 



and an operator G that causes the inversion of sign of 6 
in Hamiltonian: 

GHe{6) = He{-5)G, [Tm, G] = 0, (B3) 

an eigenvector |^) of He with the eigenvalue E{S^ F) sat- 
isfies: 

He{T2i\^)} = {E{5,F) - 2lF}{T2i\^)} (B4) 
He{T2i+i\^)} = {E(-5,F) - {21 + l)F}{T2i+i\^)}. 

Thus, the eigenenergies of Hamiltonian 



Eo,„ = E{S,F) + 2nF 

Ei,n = E{-5,F) + {2n+l)F, 



(B5) 



consists of two Wannier- Stark ladders with the corre- 
sponding eigenstates \^f3,n) = T_(2n+/3)G'^|^) [22,]. A 

I 



further calculation can prove that Eq = E{S^ F) = 
-E{-5,F). 

For an initial state expanded in Wannier- Stark basis: 

1^) = ^Co,n|^0,n) + (B6) 



the dynamics of |$) under Hamiltonian is given by 



\m) = E co,„e-'^«."*|*o,„> + E ci,ne-'^^-*|*i,n). 



(B7) 



Expanding Wannier- Stark functions in Bloch basis: 



,n/ — / (^iSn 



(B8) 



and projecting \^{t)) onto Bloch basis, one can get 



{Xo{k)mt)) = e-^^»*[ao,o(fc)Co(fc + Ft) + ai,o(fc)e-^(^-2^»)*Ci(fc + Ft)], 
{Xi{k)\m) = e-*^°*[6o,o(fc)Co(fc + Ft) + h,o{k)e-'^^-^^°'>'Ci{k + Ft)], 

I 



(B9) 



where C(3 are the Fourier series of Cfs^n- 



riod of 



(BIO) 



To = 



2n 



F-2Eo 



(B12) 



which are 7r-periodic^ To get Eq.(jB9l) one has to Ti is half of the Bloch time Tb = 2%/F for the single 

use T-2n\xp{k)) = e '^"^[^^(fc)) (translation of Bloch band system((5 = 0). In general if Ti and are com- 

waves). From Ea. (|B9|l . one can see that the dynamics of mensurate 
a particle is characterized by two periods: Cp are func- 
tions with period of 



(Bll) 



Ti 2F m 

T2 = F-2Eo=n ^^^^^ 



thus the wavefunction is reconstructed at integer multi- 
whereas the exponential function g-*(^-2^o)t j^g^g ^ pg_ p}gg of Bloch-Zener time (Tbz = nTi). 
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